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ABSTRACT 

The  aim  of  this  paper  is  to  simulate  electro-thermal  effects  in  Heterojunction  Bipolar  Transistors  (HBTs)  based 
on  Si/SiGe  and  on  AlGaAs/GaAs,  by  means  of  an  analytical  electro -thermal  model,  already  proposed  by  us, 
able  to  calculate  the  temperature  and  current  distribution  for  any  integrated  device,  whose  structure  can  be 
represented  as  an  arbitrary  number  of  superimposed  layers  with  a 2-D  embedded  thermal  source. 
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I.  Introduction 

One  of  the  main  differences  between  Si  based  VLSI  and  RF  electronics  is  the  choice  of 
semiconductor  materials  and  transistor  types.  While  Si  is  the  only  semiconductor  used  in  VLSI,  a 
wide  range  of  alternative  materials  and  devices  are  present  in  RF  electronics.  The  general  evolution  of 
these  electronic  devices  emphasizes  the  growing  importance  of  the  thermal  problem  during  the  design 
process,  where  analytical  models  for  the  temperature  evaluation  are  useful  tools  to  calculate  the 
optimal  set  of  geometrical  parameters  that  minimize  the  device  thermal  phenomena. 

Heterojunction  bipolar  transistors  (HBTs)  have  attained  enough  maturity  as  RF  power  devices  due  to 
their  intrinsic  high-power  density,  linearity  and  efficiency  [1],  Moreover  HBTs  are  actually  used  in 
very  high  speed  applications.  They  are  capable  of  microwave  performance  up  to  100  GHz.  In  order  to 
push  the  frequency  performance  up  to  its  limit,  the  operating  point  of  the  HBT  lies  at  relatively  high 
current  density,  and  this  current  density  increases  with  used  semiconductor  technologies.  Therefore  it 
is  important  to  study  the  electro-thermal  performance,  because  the  self  heating  appears,  mainly  in  the 
base-collector  depletion  region,  where  both  current  and  field  are  high. 

The  analytical  electro-thermal  modelling  of  electronic  devices  is  such  a difficult  problem  to  deal  with 
that  the  analytical  models  that  have  been  proposed  until  now  are  either  over-simplified  or  rather 
inefficient  from  a computational  point  of  view.  The  reason  for  that  lies  in  the  complex  structure  of  an 
integrated  device  and  in  the  non-linear  thermal  properties  of  the  materials.  However,  only  a physical- 
based  analytical  model  can  give  the  proper  physical  insight  in  order  to  understand  the  connections 
between  a number  of  geometrical  and  technological  parameters  and  the  device  electro -thermal 
performance.  Unfortunately,  an  analytical  model  suffers  the  unavoidable  simplifying  hypothesis  by 
which  the  numerical  calculation  can  be  carried  out.  In  spite  of  that,  the  analytical  model  for  the 
temperature  evaluation  is  a useful  tool  during  the  design  process  to  calculate  the  optimal  set  of 
geometrical  parameters  that  minimise  the  thermal  phenomena  in  an  integrated  device.  Furthermore,  in 
case  of  multifinger  devices  or  thermal  coupling  between  contiguous  devices,  the  inaccuracy 
introduced  by  the  ipothesis  of  uniform  channel  temperature  can  be  relevant. 
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In  this  paper  we  present  a simulation  study  of  electro-thermal  effects  in  HBTs  based  on  Si/SiGe  and 
on  AlGaAs/GaAs  in  order  to  compare  them,  obtained  by  an  analytical  electro-thermal  model,  already 
proposed  by  us  [2-5],  In  particular  the  problem  of  the  self-heating  of  transistors  fabricated  on  GaAs 
substrates  is  investigated  since  GaAs  has  a poor  thermal  conductivity.  The  electro-thermal  feedback  is 
considered  evaluating  the  output  current  as  a function  of  the  local  temperature  in  the  active  region. 

The  calculation  times  for  the  proposed  simulations  are  quantified  in  a minute  for  the  main  part  of  the 
algorithm,  and  in  less  than  0.01s  for  the  graphic  routines. 

The  presentation  of  the  paper  is  organized  as  follows.  In  Section  II  we  briefly  review  our  electro- 
thermal model,  based  on  the  solution  of  the  3-D  steady-state  heat  equation  with  temperature- 
dependent  thermal  conductivity  for  a single  integrated  device  or  a given  configuration  of  two  or  more 
devices.  Then,  in  Section  III  two  application  examples  of  Si/SiGe  and  AlGaAs/GaAs  HBTs  are 
given,  together  with  the  discussion  of  relative  results.  The  conclusions  are  described  in  Section  IV. 

II.  A Brief  Review  of  our  Electro-Thermal  Model 

An  exhaustive  description  of  our  model  is  in  [2-5].  In  this  Section  we  briefly  describe  our  model,  able 
to  calculate  the  temperature  and  current  distribution  for  any  integrated  device,  whose  structure,  shown 
in  Fig.  1,  can  be  represented  as  an  arbitrary  number  of  superimposed  m layers  with  a 2-D  embedded 
thermal  source,  so  as  to  include  the  effect  of  the  package. 


Figure  1.  Typical  electron  device  including  coating,  die  attachment,  mounting  and  heat  sink  (from  [2]). 

For  a complete  model,  the  contribution  of  the  mounting  and  die-attachment  layers  (see  layers  (c)  and 
(d)  in  Fig.  1 to  the  temperature  rise  and  thermal  resistance  has  been  taken  into  account.  The  top  path 
for  the  heat  flux  has  also  been  included,  e.g.  mould  compound  (see  layer  (g)),  where  natural 
convection  between  the  package  top  surface  and  the  still  air  is  the  mechanism  of  heat  loss. 

The  aims  of  our  model  have  been: 

• to  solve  analytically  the  non-linear  3-D  steady-state  heat  equation 

• to  take  into  account  the  dependence  on  temperature  of  the  thermal  conductivity 

• to  take  into  account  the  interaction  between  two  neighbouring  devices 

• to  implement  the  coupling  between  the  electrical  and  thermal  behaviour  of  the  device  to 
determine  the  actual  channel  temperature 

• to  take  into  account  the  presence  of  the  package 

• to  develop  a full  general  model  useful  to  simulate  any  multilayer  electronic  and  optoelectronic 
device 
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• fastness  and  accuracy 

• to  implement  it  on  a Personal  Computer. 

In  particular,  for  the  structure  of  Fig.  1,  in  order  to  determine  the  temperature  distribution,  the 
following  non-linear  steady-state  heat  equation  has  to  be  solved  [1]: 

V • [kTH(T)VT(x,y,z)]  = -Qv(x,y,z)  (1) 


where  T(x,y,z ) is  the  temperature  field,  km  is  the  temperature -dependent  thermal  conductivity  and 
Qv(x,y,z ) is  the  dissipated  power  density. 

The  basic  assumptions  of  the  model  have  been: 

1.  the  device  and  package  structures  can  be  represented  as  a set  of  superimposed  homogeneous 
layers; 

2.  the  thickness  of  each  layer  is  constant; 

3.  the  extension  of  the  layers  in  the  x and  y directions  is  infinite; 

4.  the  contact  thermal  resistance  is  neglected; 

5.  the  thermal  source  is  modeled  as  a 2-D  geometrical  shape  Qs(x,y),  located  at  the  interface  between 
two  contiguous  layers,  say  the  k- th  and  the  (k+\  )-th; 

6.  the  device  self-heating  is  due  to  the  Joule  heating  and  other  contributions  are  neglected. 

Eq.  1 can  be  solved  considering  the  following  expression: 

V-[kXH(T)VT(x,y,z)]  = 0 (2) 

and  accounting  for  the  heat  source  in  the  Boundary  Conditions  (BCs). 

We  have  applied  the  Kirchhoff  transformation  to  each  layer  in  the  following  form: 


A6(x,y,z) 


T(x,y,z)-To 


J k(x)dx 


T0 


(3) 


so  as  to  linearise  Eq.  (2): 


V20(x,y,z)  = 


-Q(x,y,z) 

k0 


(4) 


where  Ad  ( x,y,z ) = 0 ( x,y,z)-To , 0 ( x,y,z ) is  the  transformed  temperature  field,  To  is  the  reference 
temperature  and  ko  is  the  constant  thermal  conductivity,  evaluated  at  a specific  reference  temperature 
To.  In  this  way  the  solution  of  the  transformed  problem  with  the  unknown  At?  is  equivalent  to  the 
solution  of  the  original  non-linear  problem  with  the  unknown  AT.  when  the  inverse  transform  is 
carried  out. 

Furthermore,  to  simplify  the  mathematical  steps,  the  2-D  Fourier  transform  has  been  applied  to  Eqn. 
(4),  leading  to  the  one -dimensional  ordinary  differential  equation  for  each  /-th  layer. 

Fig.  2 shows  the  flow-chart  of  the  main  mathematical  steps  involved  in  our  method  [2-5]. 

The  electro-thermal  feedback  is  considered  by  evaluating  the  output  current  as  a function  of  the  local 
temperature  in  the  active  region  and  iterating  the  solution  of  (4)  using  the  actual  value  of  Q(  6). 


270 


Vol.  9,  Issue  3,  pp.  268-277 


International  Journal  of  Advances  in  Engineering  & Technology,  June,  2016. 

©IJAET  ISSN:  22311963 


Figure  2.  Flow-chart  of  the  mathematical  steps  involved  in  the  proposed  method  (from  [3]). 


2.a  Electro-thermal  model  applied  to  HBT 

In  order  to  solve  Eq.  (4),  we  are  referred  to  a typical  HBT  structure  [1],  in  which  we  have  used  the 
name  AlGaAs/GaAs  HBT  for  a device  with  GaAs  as  the  base  material  and  AlGaAs  as  the  emitter. 
This  is  also  the  same  for  SiGe  based  transistors. 

In  particular  we  have  represented  the  device  as  a surface  without  any  thickness,  where  X and  Y are 
the  dimensions  of  the  periodic  cell,  in  which  we  have  divided  the  device. 

The  overall  heat  source,  which  is  located  in  the  depletion  layer  at  the  collector,  has  been  divided  into  a 
set  of  elementary  point  sources  located  in  the  middle  of  each  elementary  XxY  periodic  cell. 

The  first  step  is  to  solve  the  electro-thermal  problem,  assuming  the  heat  generation  and  the  current  as 
mutually  dependent  for  each  single  elementary  device  which  is  related  to  a single  elementary  heat 
source;  finally  the  resulting  thermal  field  of  the  whole  structure  is  obtained  by  the  superposition  of  all 
the  elementary  fields. 

The  solution  of  (4),  for  a point  thermal  source  can  be  expressed  as  [6]; 


A9(x,  y,z)  = 


Q(x0i , Yoi  ’ z0i ) 


27Ut0  V(x  - x0i  ) + (y  - Yoi ) + (z  - z0i  Y 


(5) 


where  xoi,  yoi  and  z.o,  are  the  coordinates  of  the  i-th  heat  point  source,  x,y,z  is  the  generic  position  in 
which  the  temperature  increase  A 0(x,y,z)  above  the  reference  is  evaluated  and  Q(xoi,yoi,Zod  is  the 
dissipated  power  of  the  i-th  elementary  heat  source. 

The  electro-thermal  feedback  can  be  implemented  for  each  XxY  subsection  considering  that  Q can  be 
assumed  as  the  temperature -dependent  electrical  power  P(6i  (xoi,yo;,Zo;))  of  the  i-th  elementary  device 
centred  in  xo;,  yo„  Zoi  and,  thus,  corresponding  to  the  i-th  heat  source  centred  in  the  same  point. 

The  electrical  power  of  the  i-th  elementary  device  is: 
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P(X0i,y0i,Z0i)  = P(0i)  - IB(0i)  Vbe  + Ic  (0i)  V CE  (6) 

where  Vbe,  the  voltage  drop  between  base  and  emitter,  is  temperature -independent,  Ib  and  7care  the 
base  and  collector  currents,  respectively,  of  the  elementary  cell  centred  in  xo„yo;,Zo;  having  temperature 
6i,  and  Vce,  the  voltage  drop  between  collector  and  emitter,  is  temperature -independent. 

The  dependence  of  the  current-voltage  equation,  based  on  physical  parameters,  on  the  temperature 
was  widely  studied  in  the  past.  The  parameters  that  have  been  taken  into  account  for  their  thermal 
dependence  are  the  electron  mobility,  saturation  velocity,  permittivity,  energy  band  gap,  threshold 
voltage  and  built-in  voltage.  The  reader  can  refer  to  for  the  empirical  expression  of  the  foregoing 
parameters  [6]. 

III.  Electro-Thermal  Model  Implementation  for  HBTs  based  on 
Si/SiGe  and  AlGaAs/GaAs:  Simulation  Results  and  Discussion 

We  have  firstly  extracted,  by  a MATLAB  code  written  by  us  only  for  academic  non-commercial 
purpose,  the  Si/SiGe  and  AlGaAs/GaAs  HBTs  1-V  curves  versus  bias  using  the  Gummel-Poon 
model,  whose  parameters  are  those  of  SPICE  simulator  [7-8]. 

3.a  HBT  based  on  Si/SiGe 

We  have  considered  a HBT  structure  having  X=  0.25  pm  and  Y=  30  pm  for  the  emitter  and  X=  50 
pm  and  Y=  300  pm  for  the  collector. 

In  Fig.  3 we  have  reported  the  Ic-Vce  characteristics  at  room  temperature  T = 300  K,  for  a base 
current  varying  from  1 mA  to  3 mA. 
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V ce(V) 


Figure  3.  Simulated  Ic-Vce  characteristics  at  T=300  K for  a base  cun'ent  varying  from  1 mA  to  3 mA. 

The  obtained  heat  generated  at  the  base-collector  junction,  by  Eq.  (6),  has  been  equal  to  0.0702  W. 
Now  it  has  been  possible  determined  the  temperature  increase  due  to  self  heating  effect. 

Fig.  4 shows  the  simulated  average  temperature  Tm  at  the  BE  junction  versus  time,  in  the  time  range 
10"83-10"3  s. 
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Figure  4.  Simulated  average  temperature  Tm  at  BE  junction  versus  time  for  Si/SiGe  HBT. 

In  the  interval  10  6-M0‘5  s,  it  clearly  appears  that  the  temperature  increases  more  rapidly,  tending  to  a 
saturation  value  equal  to  7.85  K. 

The  electro-thermal  feedback  is  shown  in  Fig.  5,  in  which  we  have  reported  the  simulated  output  I-V 
characteristics  both  at  300  K (blue  lines)  and  at  307.85  K (red  lines).  This  figure  shows  the  well- 
known  degradation  effect  on  I-V  curves  due  to  temperature  increase. 
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Figure  5.  Simulated  I-V  characteristics  both  at  300  K (blue  lines)  and  at  307.85  K (red  lines). 

We  have  also  considered  a Si/SiGe  HBT,  whose  structure  is  the  same  of  the  previous  one,  except  for 
the  dimension  X of  the  emitter,  equal,  in  this  case,  to  2.5  pm. 

Fig.  6 shows  the  simulated  average  temperature  Tm  at  the  base-emitter  (BE)  junction  versus  time.  For 
this  device  the  saturation  value  of  Tm  is  4.55  K and  this  result  allows  to  say  that  the  self  heating  is  in 
inverse  proportion  to  the  emitter  area. 

Moreover,  as  shown  in  Fig.  7,  in  which  we  have  reported  the  simulated  I-V  curves  both  at  300  K (blue 
lines)  and  at  304.55  K (red  lines),  the  device  feels  more  the  effects  of  temperature  increase  at  high 
base  currents,  with  a resulting  decrease  of  device  efficiency. 
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Figure  6.  Simulated  average  temperature  Tm  at  BE  junction  versus  time  (X  = 2.5  pm  for  the  emitter). 


Figure  7.  Simulated  I-V  characteristics  both  at  300  K (blue  lines)  and  at  304.55  K (red  lines). 

3.b  HBT  based  on  AlGaAs/GaAs 

We  have  repeated  the  same  simulations,  considering  now  an  AlGaAs/GaAs  HBT,  having  X=  0.25  pm 
and  Y=  30  pm  for  the  emitter  and  X=  50  pm  and  Y=  300  pm  for  the  collector. 

In  Fig.  8 we  have  reported  the  Ic-Vce  characteristics  at  room  temperature  T = 300  K,  for  a base 
current  varying  from  1 mA  to  3 mA.  The  obtained  heat  generated  at  the  base-collector  junction, 
calculated  by  Eq.  (6),  has  been,  for  this  device,  equal  to  0.0508  W. 
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Figure  8.  Simulated  Ic-Vce  characteristics  at  T=300  K for  a HBT  based  on  AlGaAs/GaAs. 
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The  temperature  Tm  versus  time  is  shown  in  Fig.  9,  in  which  it  clearly  appears  that  Tm  tends  to  a 
saturation  value  equal  to  8.4  K. 


time(s) 


Figure  9.  Tm  at  BE  junction  versus  time  for  AlGaAs/GaAs  HBT. 


The  electro-thermal  feedback  is  shown  in  Fig.  10,  in  which  we  have  reported  the  simulated  output  I-V 
characteristics  both  at  300  K (blue  lines)  and  at  308.4  K (red  lines). 


Figure  10.  Simulated  I-V  characteristics  both  at  300  K (blue  lines)  and  at  308.4  K (red  lines). 

Also  for  this  device,  we  have  modified  the  emitter  dimension  (X  = 2.5  pm),  obtaining  that  the 
saturation  value  of  Tm  becomes  5.4  K,  as  shown  in  Fig.  1 1 . 


L 

r 

j 

io's  io'7  io'6  io'5  io’4  10  s 

time(s) 


Figure  11.  Tm  at  BE  junction  versus  time  for  AlGaAs/GaAs  HBT  (X  = 2.5  pm  for  the  emitter). 
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In  Fig.  12  we  have  reported  the  temperature  Tm  versus  time  for  two  HBTs  having  the  same 
dimensions,  but  different  substrates  (Si,  red  line  and  GaAs,  blue  line). 
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Figure  12.  Tm  versus  time  for  two  HBTs  having  different  substrates  (Si,  red  line  and  GaAs,  blue  line). 

It  is  easy  to  observe  that  Si  substrate  allows  to  have  smaller  values  of  Tm,  with  a transitory  time 
shorter. 

Considering  the  Early  effect,  the  breakdown  voltage  values  in  Si/SiGe  HBT  are  relatively  low  due  to 
their  smaller  band-gap  of  the  collector  material.  This  makes  it  more  difficult  to  use  them  for  high- 
power  applications,  where  AlGaAs/GaAs  HBTs  are  preferred.  However,  this  problem  can  be  reduced 
by  the  DHBT  approach  since  designing  the  collector-base  junction  as  a heterojunction  allows  to  use  a 
wide  band-gap  material  as  the  collector  material  leading  to  high-breakdown  values  [9-11]. 

The  calculation  time  for  the  proposed  examples  has  been  quantified  in  a minute  for  the  main  part  of 
the  algorithm  and  in  less  than  0.01  s for  the  graphic  routines. 

All  simulations  were  carried  out  in  ADS  2014  on  an  Asus  K55VD  computer  which  uses  an  Intel  Core 
i-7  3630QM  processor  running  at  2.4  GHz,  with  4 GB  of  RAM. 

IV.  Conclusions 

In  this  paper  we  have  presented  a simulation  study  of  electro-thermal  effect  on  performance  of  HBTs 
based  on  Si/SiGe  and  AlGaAs/GaAs,  by  means  of  an  analytical  electro-thermal  model,  already 
proposed  by  us,  able  to  calculate  the  temperature  and  current  distribution  for  any  integrated  device, 
whose  structure  can  be  represented  as  an  arbitrary  number  of  superimposed  layers  with  a 2-D 
embedded  thermal  source.  The  electro-thermal  feedback  has  been  considered  by  evaluating  the  output 
current  as  a function  of  the  local  temperature  in  the  active  region. 

Moreover  it  has  been  possible  determined  the  temperature  increase  due  to  self  heating  effect,  in  order 
to  evaluate  the  average  temperature  Tm  at  the  BE  junction.  In  this  way  we  have  demonstrated  that  Si 
substrate  allows  to  have  smaller  values  of  Tm,  with  a transitory  time  shorter. 
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